Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1326 1500 (N32W05 L=123)   345     94        0 (1294/1454-1354/1539) -33.7 W01-W07 N29-N34
      2 1674 1424 (N39W35 L=153)   915    214        2 (1580/1374-1763/1473)   2.9 W25-W45 N37-N40
      3 1662 1264 (N28W35 L=153)  1234    245        2 (1577/1159-1723/1327) -19.7 W26-W41 N22-N30
      4  583 1277 (S01E33 L= 86)   709    149        0 ( 507/1226- 630/1314)   0.0 E38-E29 S03-N00
      5 1484 1209 (N19W22 L=141)   684    260        2 (1395/1131-1557/1304)  65.0 W18-W25 N11-N26
      6  362 1131 (S15E45 L= 73)   333    157        0 ( 308/1088- 427/1152)  33.7 E49-E40 S18-S12
      7 1813 1001 (N17W54 L=172)   966    131        3 (1776/ 949-1849/1072) -40.6 W48-W58 N14-N20
      8  195  979 (S27E61 L= 57)   473     53        0 ( 173/ 950- 221/1000)  29.7 E64-E57 S29-S25
      9  187  889 (S33E63 L= 56)  1596    187        5 ( 155/ 830- 235/ 944) -45.0 E69-E55 S36-S29
     10  815  892 (S17E08 L=110)   196     91        0 ( 798/ 851- 834/ 937)  90.0 E09-E07 S19-S14
     11 1679  777 (N00W48 L=166)   740    151        1 (1654/ 700-1703/ 843) -76.0 W45-W50 S04-N04
     12  173  791 (S40E68 L= 50)  1484    121        2 ( 129/ 765- 223/ 819)   8.4 E85-E58 S41-S36
     13 1512  779 (S04W36 L=154)   496     97        0 (1483/ 745-1544/ 814)  -9.5 W32-W38 S05-S02
     14 1307  743 (S12W24 L=143)   261     78        0 (1275/ 720-1339/ 773) -11.3 W21-W26 S13-S11
     15  753  703 (S30E08 L=111)   679    213        2 ( 682/ 632- 841/ 764)  59.7 E11-E03 S36-S24
     16 1358  309 (S35W53 L=171)   598     98        1 (1317/ 293-1412/ 333)  12.5 W47-W57 S36-S33
     17 1541  294 (S29W73 L=191)  1020     88        1 (1506/ 267-1580/ 327)   7.6 W64-W82 S30-S27
  

    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