Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  859 1426 (N31E09 L=121)   551    197        0 ( 802/1350- 893/1497) -63.4 E13-E07 N25-N35
      2 1253 1345 (N23W17 L=147)   578    164        0 (1188/1302-1309/1369) -14.0 W12-W20 N20-N24
      3  232 1297 (N23E63 L= 67)  1125    120        3 ( 205/1245- 261/1364)  90.0 E66-E60 N19-N27
      4 1813 1260 (N13W62 L=192)   514     77        2 (1784/1236-1847/1288) -15.9 W58-W66 N11-N14
      5  172 1244 (N19E71 L= 59)  1065    107        1 ( 150/1191- 190/1287)  63.4 E75-E67 N16-N22
      6 1536 1194 (N12W35 L=165)   314     86        3 (1513/1166-1571/1223) -45.0 W32-W37 N10-N13
      7 1926  906 (S11W76 L=206)  1215    120        2 (1912/ 851-1944/ 972)  69.4 W72-W80 S14-S06
      8 1825  930 (S07W58 L=189)   625     79        1 (1811/ 881-1846/ 968)  90.0 W56-W60 S10-S04
      9 1623  842 (S11W39 L=169)   443    112        0 (1598/ 789-1642/ 889)  71.6 W37-W40 S13-S07
     10  929  742 (S12E07 L=123)   432    211        1 ( 861/ 677- 991/ 809)  37.9 E12-E03 S16-S08
     11 1511  722 (S17W31 L=161)  1121    342        4 (1408/ 649-1592/ 776)  -8.1 W23-W37 S21-S13
     12 1797  700 (S22W59 L=189)  1322    173        1 (1738/ 654-1854/ 770)  24.8 W52-W65 S24-S18
     13 1560  656 (S22W36 L=166)   303     59        1 (1533/ 627-1587/ 676)  14.0 W33-W37 S24-S21
     14 1498  597 (S26W31 L=161)   276     72        1 (1480/ 573-1521/ 632)  90.0 W29-W33 S27-S23
     15  577  532 (S25E35 L= 95)   934    147        2 ( 555/ 464- 601/ 598)  59.0 E37-E31 S30-S20
     16  880  457 (S32E14 L=116)   247    104        0 ( 845/ 430- 904/ 501)  51.3 E16-E11 S33-S28
     17  470  434 (S33E50 L= 80)   986    159        4 ( 404/ 414- 561/ 450)  -3.6 E57-E40 S34-S31
  

    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