Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  821 1730 (N57E10 L= 75)   481     58        2 ( 800/1693- 845/1766)  81.9 E12-E06 N53-N60
      2 1779 1324 (N13W61 L=145)   975    194        1 (1723/1245-1834/1399) -61.4 W56-W65 N07-N18
      3 1436 1348 (N20W32 L=116)   297    128        0 (1368/1335-1481/1368)   0.0 W26-W35 N19-N21
      4  136 1129 (N17E75 L=  9)   881     72        1 ( 107/1098- 156/1152)  -7.6 E85-E69 N15-N18
      5 1899 1006 (S08W68 L=153)  2620    158        2 (1848/ 940-1939/1079) -55.0 W61-W78 S12-S02
      6  213  913 (N04E62 L= 23)   434     56        1 ( 192/ 885- 236/ 937)  26.6 E64-E58 N02-N05
      7 1807  896 (S13W56 L=140)   600    138        2 (1737/ 870-1843/ 934)   0.0 W48-W60 S14-S10
      8 1348  725 (S17W18 L=102)   487    123        1 (1319/ 662-1383/ 778)  60.3 W15-W20 S20-S13
      9  963  542 (S24E09 L= 75)   248     67        0 ( 942/ 517- 983/ 564)  45.0 E11-E07 S26-S22
  

    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