Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  784 1413 (N25E11 L= 50)   847    202        1 ( 675/1364- 855/1441) -22.6 E18-E06 N21-N28
      2 1139 1320 (N16W10 L= 71)   265    126        0 (1094/1303-1206/1345) -23.2 W07-W14 N14-N17
      3  231 1250 (N22E56 L=  5)  1684    287        1 ( 152/1163- 344/1337) -24.9 E71-E43 N14-N27
      4  652 1022 (N03E22 L= 39)   266     62        1 ( 616/ 999- 683/1040) -26.6 E24-E20 N01-N04
      5 1726  915 (S14W46 L=107)   307     52        3 (1703/ 890-1754/ 937)  14.0 W44-W48 S15-S13
      6 1660  734 (S25W41 L=103)   923    223        0 (1635/ 638-1683/ 855)  85.9 W38-W44 S31-S17
      7 1753  664 (S31W52 L=114)   826    107        1 (1706/ 638-1792/ 707)  14.0 W47-W57 S32-S27
      8  364  493 (S26E59 L=  2)  1843    225        3 ( 324/ 381- 402/ 585)  30.6 E75-E51 S32-S19
  

    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