Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  439 1589 (N20E58 L=257)   587     88        1 ( 409/1540- 467/1624) -80.5 E59-E55 N16-N22
      2  658 1571 (N23E38 L=276)   316     80        0 ( 648/1528- 672/1608) -51.3 E40-E36 N20-N25
      3 1007 1514 (N27E13 L=302)  1107    342        1 ( 864/1416-1128/1576) -10.3 E24-E02 N23-N29
      4 1209 1473 (N29W01 L=316)   230     86        0 (1179/1446-1244/1510) -31.0 E01-W04 N28-N31
      5  555 1412 (N11E39 L=276)   302     61        2 ( 529/1383- 577/1443) -63.4 E40-E36 N09-N13
      6 1302 1333 (N22W10 L=325)   223    105        0 (1269/1300-1338/1373) -26.6 W07-W13 N21-N24
      7 1762 1227 (N27W49 L=  3)  1699    210        3 (1729/1127-1789/1330)  70.3 W45-W53 N20-N34
      8  287 1129 (S11E53 L=262)   631     85        1 ( 254/1096- 321/1155)  38.7 E56-E49 S12-S08
      9 1597 1048 (N12W35 L=350)   260     46        1 (1579/1022-1622/1071) -45.0 W33-W37 N10-N12
     10 1877  952 (N13W65 L= 20)  1973    167        5 (1832/ 878-1922/1007)  66.0 W60-W74 N07-N16
     11  361  925 (S22E43 L=271)   938    178        1 ( 295/ 887- 444/ 960)   8.1 E49-E35 S24-S19
     12 1809  915 (N09W57 L= 11)   501     66        2 (1794/ 876-1826/ 950) -68.2 W55-W59 N06-N11
     13  415  844 (S25E37 L=277)   308     91        3 ( 394/ 799- 444/ 880)  63.4 E39-E34 S28-S22
     14 1250  620 (S21W24 L=338)  1362    342        1 (1195/ 494-1354/ 712)  78.7 W21-W28 S28-S13
     15  335  605 (S44E49 L=266)   666     83        1 ( 316/ 564- 352/ 645)  49.4 E51-E45 S47-S40
     16  726  591 (S35E10 L=304)   354    170        0 ( 658/ 570- 801/ 624)  26.6 E15-E05 S37-S31
     17 1368  489 (S26W37 L=351)   628    194        1 (1271/ 472-1451/ 504)  17.1 W29-W42 S28-S23
     18 1526  402 (S27W56 L= 10)   414     62        0 (1511/ 374-1540/ 442) -45.0 W52-W58 S29-S24
  

    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