Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1102 1388 (N27W07 L=334)  1230    346        1 ( 922/1331-1250/1433) -20.9 E05-W16 N22-N30
      2 1565 1276 (N17W38 L=  6)   367     59        1 (1551/1235-1585/1316)  59.0 W36-W39 N14-N19
      3 1386 1260 (N17W25 L=352)  1013    198        4 (1321/1198-1454/1313) -53.1 W20-W29 N12-N20
      4  658 1177 (N15E23 L=304)   510    181        0 ( 621/1089- 692/1264) -74.7 E25-E20 N09-N20
      5 1779  875 (S10W53 L= 21)  1516    303        1 (1702/ 782-1845/ 952) -40.6 W46-W61 S16-S04
      6  425  904 (S01E41 L=287)   346     91        2 ( 387/ 883- 456/ 924)   0.0 E43-E37 S02-N00
      7 1399  813 (S11W23 L=350)   386     82        0 (1377/ 766-1423/ 860)  71.6 W21-W24 S13-S07
      8 1226  639 (S21W11 L=339)   649    216        0 (1132/ 574-1333/ 729)  34.7 W05-W18 S24-S15
      9 1490  564 (S28W31 L=358)  2262    451        4 (1308/ 439-1642/ 722)  39.3 W18-W42 S36-S18
  

    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