Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  348 1441 (N13E56 L=160)   556     87        2 ( 320/1394- 374/1472)  90.0 E58-E54 N09-N15
      2 1785 1243 (N26W51 L=267)   399     57        1 (1761/1218-1808/1265) -45.0 W48-W53 N24-N27
      3  414 1184 (S01E42 L=174)   305     94        0 ( 380/1169- 454/1202)   0.0 E45-E39 S02-N00
      4 1901 1078 (N18W66 L=282)  1067    139        1 (1868/1033-1928/1111) -26.6 W61-W70 N15-N19
      5 1311  908 (S03W19 L=235)   311     85        0 (1280/ 870-1340/ 942) -31.0 W16-W21 S04-N00
      6 1324  773 (S11W23 L=239)   237     74        2 (1310/ 734-1337/ 819) -78.7 W21-W23 S12-S07
      7  894  542 (S33W00 L=217)   365    109        0 ( 838/ 522- 953/ 570)   0.0 E03-W04 S33-S31
      8 1293  537 (S26W27 L=244)   571     88        0 (1246/ 511-1335/ 561)  -7.1 W23-W31 S26-S24
      9 1736  503 (S19W66 L=282)   929    250        1 (1699/ 483-1776/ 526)   0.0 W60-W73 S20-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