Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1215 1830 (N61W38 L=211)   505     64        0 (1201/1790-1230/1862)  26.6 W32-W46 N56-N64
      2 1818 1319 (N15W66 L=239)   456     51        0 (1788/1299-1845/1334) -18.4 W61-W69 N13-N15
      3 1340 1274 (N18W23 L=196)   405    187        0 (1279/1238-1410/1332) -36.9 W19-W27 N15-N21
      4  526 1092 (N12E32 L=141)   724    227        1 ( 498/1014- 570/1192) -65.6 E34-E29 N07-N18
      5 1803  978 (S05W56 L=229)  1083    127        1 (1749/ 937-1860/1009) -14.0 W50-W62 S07-S02
      6  533  959 (N04E32 L=141)   296     60        0 ( 511/ 927- 553/ 987)  53.1 E33-E30 N01-N05
      7  524  694 (S13E36 L=138)   601    159        2 ( 474/ 665- 582/ 735)  16.7 E40-E30 S14-S10
      8 1161  555 (S26W06 L=179)  1896    523        4 (1012/ 377-1334/ 729)  45.0 E06-W17 S38-S15
      9 1635  606 (S28W42 L=216)   596    123        0 (1600/ 556-1675/ 659)  50.2 W39-W45 S30-S24
  

    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