Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1183 1748 (N54W04 L= 84)   404     52        2 (1153/1727-1218/1764)   0.0 W01-W08 N52-N55
      2  465 1522 (N27E49 L= 30)  1350    227        0 ( 419/1432- 497/1633) -63.4 E56-E46 N20-N34
      3 1351 1390 (N28W19 L= 98)  1525    462        0 (1184/1250-1436/1530) -40.1 W05-W24 N19-N35
      4 1778 1393 (N32W60 L=140)   485     70        1 (1752/1362-1798/1421) -63.4 W57-W63 N29-N33
      5 1138 1228 (N15W05 L= 85)   394    168        0 (1103/1179-1193/1294) -39.8 W02-W08 N12-N18
      6 1575 1208 (N18W35 L=115)   380    153        0 (1544/1177-1602/1252) -76.0 W33-W37 N15-N20
      7  440  727 (S23E38 L= 41)   283     57        0 ( 408/ 704- 468/ 751)  26.6 E41-E35 S25-S21
      8 1256  679 (S18W19 L= 98)   573    173        0 (1187/ 629-1327/ 730)  48.8 W14-W22 S21-S13
      9 1684  603 (S19W55 L=134)  1250    202        2 (1652/ 515-1727/ 674) -79.7 W51-W59 S24-S13
     10 1616  572 (S21W49 L=129)   390     51        1 (1593/ 551-1640/ 603)  71.6 W47-W51 S22-S19
     11  422  509 (S39E46 L= 34)  1081    173        1 ( 371/ 475- 499/ 553)  24.0 E53-E35 S42-S34
  

    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