Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  532 1616 (N23E51 L= 41)   628    134        2 ( 485/1594- 596/1651)  -4.1 E57-E43 N21-N24
      2 1056 1313 (N16E04 L= 88)   808    342        8 ( 978/1224-1147/1381) -18.4 E10-W02 N11-N19
      3 1298 1314 (N21W10 L=103)   342     45        2 (1271/1288-1321/1343)  76.0 W08-W11 N19-N23
      4 1304 1235 (N16W12 L=105)   609    190        2 (1248/1168-1348/1266) -26.6 W08-W14 N12-N18
      5  275 1006 (S18E52 L= 40)   378     77        1 ( 254/ 970- 308/1038) -53.1 E54-E48 S20-S16
      6  949  604 (S29W05 L= 97)  1231    437        6 ( 764/ 539-1125/ 702)  43.5 E06-W13 S37-S18
      7 1701  649 (S09W56 L=148)   365     40        2 (1678/ 627-1728/ 673)  26.6 W53-W58 S10-S07
      8  812  363 (S50W02 L= 95)   373     58        1 ( 789/ 338- 840/ 387) -11.3 W00-W05 S52-S48
      9 1091  299 (S47W29 L=121)  1630    269        2 ( 981/ 272-1207/ 351)  26.6 W19-W38 S50-S41
  

    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