Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1407 1715 (N46W22 L=260)  1608    211        1 (1315/1655-1498/1753) -10.6 W13-W30 N42-N49
      2  300 1518 (N19E63 L=175)  2416    187        6 ( 263/1432- 332/1604) -56.3 E69-E57 N13-N25
      3 1243 1408 (N20W09 L=247)  1204    300        1 (1122/1310-1352/1502) -29.4 W00-W16 N15-N25
      4 1824 1301 (N23W56 L=294)  2346    430        2 (1733/1120-1883/1445) -62.2 W48-W66 N12-N31
      5 1662 1303 (N20W39 L=277)   749    121        1 (1620/1261-1714/1344) -20.6 W36-W44 N17-N22
      6 1759 1310 (N22W49 L=287)   785    110        3 (1707/1280-1802/1344) -14.0 W43-W53 N20-N23
      7  150 1031 (S12E67 L=171)   718     80        3 ( 127/ 999- 178/1074) -45.0 E70-E62 S13-S09
      8  199 1039 (S11E60 L=178)   374     52        1 ( 181/1009- 213/1067)  90.0 E62-E58 S13-S09
      9  344 1018 (S12E45 L=192)   398     92        1 ( 314/ 975- 378/1066)  51.3 E47-E42 S14-S09
     10  746  879 (S17E16 L=222)   337    169        1 ( 676/ 852- 824/ 904)  16.7 E20-E10 S19-S15
     11  230  845 (S23E58 L=180)   461     72        1 ( 192/ 828- 261/ 874)  14.0 E62-E54 S23-S21
     12  516  706 (S30E32 L=206)   628    171        0 ( 463/ 649- 592/ 770)  36.9 E37-E25 S34-S25
     13 1558  721 (S16W39 L=277)   283     32        1 (1540/ 700-1579/ 750) -33.7 W37-W40 S17-S14
     14 1374  573 (S28W30 L=268)  1902    428        0 (1224/ 462-1539/ 686)  51.7 W22-W39 S37-S18
     15 1705  593 (S21W57 L=295)   593    112        0 (1662/ 547-1744/ 639)  63.4 W54-W59 S24-S17
     16  401  520 (S44E48 L=189)   346     50        0 ( 385/ 488- 419/ 541)  39.8 E51-E45 S46-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