Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1037 1918 (N62W62 L=235)   609     47        1 (1013/1890-1055/1943)   5.2 W53-W75 N60-N63
      2  940 1844 (N61W37 L=210)   741     97        2 ( 891/1813- 983/1876)  -3.8 W29-W46 N58-N63
      3  779 1782 (N64W13 L=186)   498     74        1 ( 748/1746- 814/1804)  21.8 W08-W18 N60-N65
      4  589 1649 (N59E16 L=157)   652     80        1 ( 566/1607- 604/1694)  76.0 E18-E13 N55-N63
      5 1593  870 (S18W30 L=203)   268     50        1 (1571/ 847-1614/ 893)   0.0 W28-W32 S19-S16
      6 1433  684 (S24W15 L=188)   485    250        0 (1376/ 563-1460/ 764)  47.7 W08-W18 S30-S19
      7  724  542 (S13E31 L=142)  1143    271        4 ( 612/ 490- 849/ 599)  -3.2 E40-E22 S14-S10
      8 1018  532 (S21E14 L=159)   247    126        0 ( 987/ 488-1054/ 591)  35.5 E17-E10 S23-S18
      9 1261  417 (S37E03 L=170)   476    111        1 (1211/ 399-1320/ 446)   0.0 E06-W01 S37-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