Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  779 1765 (N40E40 L=305)   436     59        2 ( 752/1736- 805/1794) -14.0 E44-E36 N38-N42
      2 1152 1594 (N36E04 L=342)   795    156        2 (1081/1544-1217/1618) -12.5 E09-W00 N34-N38
      3  782 1426 (N17E23 L=323)   214     99        0 ( 749/1396- 810/1472)  68.2 E24-E22 N15-N20
      4  137 1118 (S12E71 L=275)   581     56        1 ( 114/1100- 163/1137)  -6.3 E75-E66 S12-S11
      5  664 1098 (S04E23 L=323)   243     96        1 ( 635/1063- 708/1132) -26.6 E25-E19 S05-S02
      6 1265  997 (N01W14 L=  0)   400    108        0 (1205/ 973-1309/1021)   0.0 W10-W17 N00-N02
      7 1731  773 (S02W53 L= 38)  1560    285        3 (1649/ 661-1801/ 893)  71.6 W48-W56 S09-N06
      8  775  752 (S22E10 L=336)   951    369        1 ( 634/ 658- 937/ 813)  39.8 E18-W00 S31-S15
  

    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