Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  369 1283 (S03E47 L= 75)  1429    299        1 ( 275/1158- 469/1372)  57.3 E53-E40 S11-N03
      2 1458 1211 (N13W22 L=144)  1105    318        1 (1411/1064-1485/1357) -56.3 W16-W26 N04-N20
      3  466 1274 (S02E39 L= 83)   246     51        2 ( 452/1239- 479/1308)  90.0 E40-E38 S04-N00
      4  150 1171 (S11E68 L= 54)  1151    192        6 ( 124/1101- 182/1226)  90.0 E71-E63 S15-S08
      5 1691 1035 (N09W41 L=163)   556    141        0 (1658/ 959-1715/1096)  83.7 W40-W42 N03-N12
      6 1860  979 (N10W58 L=180)  1614    230        5 (1818/ 894-1893/1056) -84.3 W55-W62 N05-N15
      7 1139  624 (S28W16 L=138)   195     56        0 (1110/ 610-1174/ 638)  36.9 W14-W18 S28-S25
      8 1732  513 (S18W64 L=186)   468     73        1 (1705/ 486-1755/ 547)  78.7 W61-W65 S20-S15
      9 1207  498 (S34W25 L=147)   325    112        0 (1156/ 483-1270/ 522)   7.1 W21-W29 S34-S31
  

    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