Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1002 1572 (N34W06 L=193)   296     63        0 ( 972/1544-1022/1594) -76.0 W04-W07 N31-N35
      2 1374 1344 (N14W26 L=212)   254    161        1 (1328/1288-1427/1377)  23.2 W22-W29 N11-N15
      3 1768 1268 (N05W54 L=241)   449     86        0 (1743/1222-1797/1315) -71.6 W52-W56 N02-N08
      4  176  993 (N08E61 L=125)   732    155        1 ( 152/ 940- 212/1057) -63.4 E63-E57 N04-N12
      5  278  921 (N03E52 L=135)   637    122        2 ( 258/ 865- 295/ 986)  50.2 E54-E49 N00-N06
      6 1303  913 (S10W15 L=202)   198     77        1 (1259/ 887-1328/ 931) -26.6 W12-W16 S11-S08
      7 1092  770 (S16W01 L=187)   783    332        2 ( 978/ 655-1235/ 889)  30.5 E07-W10 S21-S10
      8  272  643 (S14E60 L=127)  1312    295        2 ( 225/ 561- 305/ 728) -79.7 E63-E55 S18-S07
      9  658  692 (S16E27 L=159)   246     65        1 ( 610/ 675- 687/ 707) -21.8 E30-E25 S16-S14
  

    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