Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1169 1617 (N33E02 L=341)   394    102        0 (1119/1603-1218/1631)  31.0 E04-W01 N30-N34
      2 1049 1496 (N22E07 L=336)   285     79        1 (1018/1470-1080/1522) -21.8 E09-E04 N20-N23
      3  983 1421 (N16E09 L=333)   219     70        0 ( 963/1403-1020/1450) -11.3 E11-E06 N15-N17
      4 1627 1235 (N17W35 L= 17)   301     85        0 (1604/1194-1653/1273) -45.0 W32-W37 N15-N19
      5 1804 1099 (N14W51 L= 34)   485     96        0 (1788/1059-1830/1144) -51.3 W49-W54 N11-N16
      6  279  960 (S21E51 L=292)  1046    213        2 ( 223/ 891- 329/1031)  48.4 E57-E46 S24-S15
      7 1170  677 (S24W16 L=358)   308    102        0 (1132/ 643-1215/ 707)  59.0 W13-W17 S26-S21
      8 1591  619 (S19W46 L= 29)   706    174        0 (1512/ 580-1653/ 663)  41.2 W41-W49 S22-S14
      9  709  582 (S39E14 L=329)   244     77        0 ( 688/ 543- 737/ 606)  68.2 E15-E11 S41-S36
     10  956  584 (S34W04 L=347)   281    108        0 ( 899/ 569-1009/ 601)  23.2 W00-W07 S35-S32
     11  244  569 (S46E69 L=274)   607     58        1 ( 226/ 545- 272/ 598) -21.8 E74-E64 S47-S43
     12  940  278 (S57W14 L=357)   607     97        1 ( 884/ 254- 980/ 308)   0.0 W07-W19 S58-S54
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y