Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1196 1565 (N35W26 L=171)   325     69        1 (1162/1537-1231/1585) -21.8 W23-W28 N33-N36
      2  556 1371 (N37E26 L=119)   900    157        2 ( 506/1314- 615/1426) -48.4 E29-E21 N32-N41
      3  321 1137 (N26E47 L= 98)  1405    343        3 ( 276/ 949- 380/1248) -73.3 E50-E43 N13-N33
      4 1697 1153 (S01W46 L=191)   368     89        1 (1679/1109-1718/1197)  36.9 W43-W48 S03-N00
      5 1759 1168 (S02W52 L=197)   503     83        0 (1725/1140-1792/1194) -38.7 W49-W54 S04-N00
      6  863 1126 (N16E08 L=137)   332    111        0 ( 818/1103- 932/1140) -15.9 E10-E03 N14-N17
      7 1074  921 (N00W01 L=145)  1007    250        1 ( 989/ 803-1122/1019) -81.3 E02-W03 S06-N07
      8  147  978 (N18E67 L= 77)   519     78        2 ( 133/ 940- 165/1016) -78.7 E68-E65 N15-N20
      9 1694  926 (S14W41 L=186)   722    117        3 (1670/ 859-1724/ 977)  60.3 W38-W43 S18-S11
     10 1039  627 (S17E08 L=137)   829    171        2 ( 964/ 552-1078/ 668) -73.3 E11-E05 S22-S12
     11  831  307 (S35E32 L=112)  1140    212        2 ( 774/ 250- 875/ 362)  26.6 E40-E26 S38-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