Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1270 1740 (N46W10 L=260)  1352    191        1 (1180/1696-1356/1773) -10.6 W01-W17 N43-N48
      2  372 1526 (N20E55 L=196)  1212    158        1 ( 355/1445- 390/1609) -47.7 E60-E50 N14-N25
      3  849 1533 (N24E18 L=233)   290     81        0 ( 816/1504- 888/1560)  51.3 E19-E15 N21-N26
      4 1057 1424 (N19E03 L=248)  1068    284        2 ( 946/1324-1163/1520) -32.7 E10-W04 N13-N24
      5 1545 1388 (N23W30 L=281)   866    266        4 (1453/1330-1656/1468) -17.4 W23-W39 N20-N28
      6 1749 1316 (N22W48 L=299)  1725    333        3 (1709/1154-1782/1456)  87.1 W45-W52 N11-N31
      7 1368 1421 (N23W18 L=268)   213     45        0 (1351/1394-1385/1449) -56.3 W16-W18 N21-N24
      8 1838 1243 (N19W57 L=308)  1552    287        1 (1798/1112-1880/1347) -77.0 W52-W62 N11-N26
      9  141  827 (S24E72 L=179)   650     52        1 ( 109/ 810- 166/ 849)   0.0 E80-E66 S24-S22
     10  346  665 (S34E49 L=201)  1251    333        0 ( 287/ 521- 432/ 786)  33.7 E65-E38 S43-S25
     11 1817  660 (S15W66 L=317)  2953    247       10 (1721/ 569-1889/ 729)   4.2 W53-W80 S20-S11
     12 1211  610 (S28W18 L=268)  1223    395        0 (1068/ 496-1368/ 713)  50.2 W10-W26 S37-S19
  

    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