Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  941 1806 (N46E34 L=348)   375    100        0 ( 912/1764- 980/1835) -76.0 E35-E32 N41-N49
      2  545 1669 (N25E54 L=328)   851     92        1 ( 502/1630- 584/1696)   0.0 E59-E48 N22-N26
      3  391 1494 (N11E56 L=326)  1125    181        3 ( 344/1429- 430/1554) -37.9 E60-E50 N07-N14
      4 1268 1388 (N24W06 L= 28)   288     80        0 (1227/1365-1309/1406)  -9.5 W03-W09 N23-N25
      5  689 1233 (N02E25 L=357)   223     78        1 ( 669/1200- 717/1267)  90.0 E26-E23 N00-N04
      6 1408 1116 (N11W20 L= 42)   279     94        0 (1385/1063-1422/1158) -59.0 W18-W22 N08-N13
      7 1841  958 (N13W59 L= 81)   877    138        0 (1796/ 889-1865/1001) -21.8 W53-W63 N09-N15
      8 1101  607 (S25W15 L= 37)   210     56        0 (1069/ 589-1128/ 626)   0.0 W12-W17 S25-S24
      9 1125  450 (S35W23 L= 45)   739    120        2 (1062/ 414-1185/ 489)   0.0 W17-W28 S35-S32
     10 1553  434 (S24W56 L= 78)   625     82        0 (1509/ 416-1596/ 452)  12.5 W50-W60 S24-S22
     11 1226  261 (S45W45 L= 67)   407     65        0 (1207/ 242-1256/ 291) -18.4 W40-W49 S46-S42
  

    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