Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  789 1792 (N34E44 L=299)   460     51        0 ( 769/1763- 811/1823) -50.2 E46-E40 N31-N37
      2 1309 1647 (N39E00 L=342)  1657    265        7 (1221/1543-1390/1731) -29.4 E09-W08 N33-N44
      3  970 1543 (N21E18 L=325)   465    144        1 ( 900/1521-1036/1558)  20.6 E21-E13 N18-N23
      4  757 1419 (N09E26 L=317)   362    137        2 ( 738/1356- 776/1480) -54.5 E28-E22 N05-N12
      5  580 1377 (N02E36 L=307)   291     61        0 ( 563/1347- 605/1410) -45.0 E37-E33 N00-N04
      6 1688 1329 (N31W35 L= 18)  1193    225        2 (1618/1259-1755/1398) -15.9 W28-W42 N26-N33
      7 1399 1166 (N11W18 L=  0)   532    181        0 (1379/1091-1438/1246)  90.0 W15-W19 N06-N16
      8 1591 1094 (N13W32 L= 15)   714    260        0 (1534/ 989-1647/1205) -34.7 W25-W38 N08-N17
      9  447 1178 (S12E40 L=303)   457    157        1 ( 378/1153- 519/1203)  21.8 E44-E34 S14-S10
     10  197 1112 (S20E63 L=280)   715    150        0 ( 148/1061- 229/1172)  37.9 E69-E59 S23-S16
     11  632  981 (S19E22 L=321)   313     86        1 ( 589/ 961- 673/1009)  31.0 E24-E19 S21-S17
     12  416  916 (S28E37 L=305)   366    135        2 ( 355/ 895- 468/ 947)  24.0 E42-E33 S30-S24
     13  574  740 (S35E22 L=321)   627    144        4 ( 508/ 711- 646/ 775)  32.5 E27-E16 S38-S31
     14 1730  543 (S11W66 L= 49)   660    104        1 (1684/ 508-1771/ 575)  50.2 W63-W69 S14-S08
     15 1341  497 (S27W37 L= 20)  1182    230        4 (1253/ 454-1439/ 569)  63.4 W32-W40 S32-S20
     16 1114  342 (S44W33 L= 16)   421     90        1 (1072/ 326-1158/ 360)  26.6 W29-W35 S45-S42
     17 1521  331 (S30W65 L= 48)   843    144        1 (1458/ 301-1588/ 357)  41.2 W60-W69 S32-S25
  

    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