Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  930 1872 (N72W05 L=152)   592     50        1 ( 912/1843- 952/1900)  45.0 W01-W10 N67-N76
      2 1596 1559 (N32W54 L=201)   554    130        0 (1549/1525-1650/1578) -20.6 W49-W59 N29-N33
      3 1123 1531 (N37W13 L=160)   354     89        0 (1074/1502-1166/1557)  14.0 W08-W16 N35-N38
      4 1198 1520 (N35W18 L=166)   264     45        1 (1182/1491-1215/1548)  45.0 W16-W20 N33-N37
      5 1556 1473 (N27W45 L=192)   696    171        1 (1520/1384-1587/1536)  78.7 W43-W47 N21-N31
      6 1703 1471 (N25W60 L=207)   674    125        3 (1661/1425-1743/1511)  18.4 W53-W65 N21-N26
      7  806 1352 (N28E12 L=135)   525    177        1 ( 770/1283- 862/1420) -68.2 E14-E08 N22-N32
      8  400 1284 (N26E44 L=103)  2520    422        7 ( 292/1172- 569/1336)  11.8 E54-E29 N18-N28
      9 1265 1264 (N18W18 L=165)   376    126        2 (1224/1225-1313/1305) -38.7 W15-W20 N15-N20
     10  399 1067 (N12E42 L=105)   624    107        2 ( 358/1025- 443/1105) -35.5 E45-E38 N09-N14
     11  501 1002 (N08E34 L=113)   277     47        1 ( 478/ 980- 530/1030)   0.0 E36-E32 N06-N09
     12 1749  858 (S13W50 L=197)   556     41        1 (1722/ 830-1771/ 885) -31.0 W47-W52 S14-S11
     13  680  796 (S06E24 L=123)   279     99        0 ( 627/ 769- 718/ 812) -18.4 E27-E21 S07-S04
     14  243  747 (S08E62 L= 85)   750     81        2 ( 222/ 697- 269/ 785)  78.7 E63-E58 S10-S05
     15  437  412 (S32E57 L= 90)   624    115        0 ( 417/ 361- 470/ 464)  63.4 E60-E53 S35-S27
     16  855  401 (S34E19 L=128)   403    126        1 ( 797/ 388- 921/ 416)  -6.3 E23-E14 S35-S32
     17 1246  249 (S54W12 L=159)   428     41        1 (1222/ 234-1273/ 271) -33.7 W08-W14 S55-S51
  

    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