Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1297 1481 (N28W29 L=184)   476    109        1 (1255/1447-1348/1515)   7.1 W25-W33 N26-N28
      2  447 1331 (N36E35 L=119)   965    253        0 ( 386/1242- 513/1412) -47.3 E41-E29 N29-N42
      3  212 1217 (N32E57 L= 97)   639     88        2 ( 188/1189- 245/1244)  -9.5 E59-E53 N30-N33
      4  629 1144 (N22E22 L=132)   219     51        0 ( 599/1127- 657/1164) -36.9 E24-E20 N20-N23
      5 1692 1123 (S03W45 L=200)   937    199        3 (1622/1098-1755/1155) -26.6 W39-W50 S05-N00
      6  231 1072 (N23E56 L= 99)  1038    207        2 ( 209/ 952- 265/1143) -85.2 E57-E51 N15-N27
      7  300  992 (N17E50 L=105)   292     82        0 ( 283/ 947- 313/1032)  51.3 E52-E48 N14-N19
      8  891  910 (N03E10 L=145)   952    201        3 ( 828/ 828- 965/1004) -79.7 E12-E06 S02-N09
      9  876  579 (S17E18 L=136)   778    196        1 ( 802/ 494- 916/ 623) -78.7 E21-E16 S22-S12
     10  714  283 (S35E42 L=112)  1106    215        3 ( 664/ 236- 760/ 329)  18.4 E49-E37 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