Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  723 1809 (N39E45 L=330)   763     85        2 ( 684/1786- 771/1835)  -9.5 E51-E39 N36-N40
      2  546 1784 (N34E60 L=314)   567     59        0 ( 533/1753- 560/1816) -36.9 E65-E55 N31-N37
      3  989 1740 (N38E19 L=355)  1139    169        1 ( 926/1674-1036/1782) -32.5 E25-E14 N33-N41
      4 1313 1671 (N41W06 L= 21)   322     74        0 (1269/1657-1352/1686)  33.7 W03-W09 N38-N42
      5  306 1549 (N15E66 L=308)   739    166        3 ( 270/1486- 342/1597) -82.9 E68-E63 N10-N18
      6  378 1458 (N09E52 L=322)   700    109        1 ( 338/1409- 417/1505)  60.3 E54-E50 N05-N12
      7  691 1439 (N12E28 L=346)   378     82        1 ( 665/1407- 724/1473) -21.8 E30-E25 N10-N14
      8  916 1406 (N14E14 L=  1)   411    124        1 ( 858/1376- 964/1426)   0.0 E17-E10 N12-N14
      9 1781 1211 (N21W48 L= 62)   428     47        0 (1748/1193-1816/1233)  14.0 W44-W52 N19-N21
     10  146 1132 (S12E69 L=305)  1155    140        2 ( 111/1080- 175/1194) -49.4 E75-E63 S15-S08
     11  207 1110 (S13E60 L=314)   921    162        2 ( 188/1027- 224/1186)  84.3 E62-E58 S18-S08
     12  186  788 (S33E65 L=309)  1367    128        2 ( 152/ 729- 214/ 852)  29.7 E75-E59 S36-S28
     13 1563  742 (S12W40 L= 54)   278    109        0 (1539/ 693-1578/ 795) -45.0 W36-W41 S14-S09
     14 1147  683 (S24W15 L= 29)   921    208        2 (1120/ 577-1179/ 788) -69.0 W12-W17 S30-S17
     15 1492  610 (S21W40 L= 54)   263     71        0 (1462/ 583-1524/ 635)  63.4 W37-W41 S22-S18
     16 1579  499 (S25W52 L= 67)  1536    305        1 (1466/ 417-1701/ 590)  56.3 W46-W58 S31-S16
     17 1332  473 (S33W35 L= 49)   705    112        2 (1273/ 447-1384/ 507)  39.8 W30-W37 S34-S29
  

    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