Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1466 1792 (N36W69 L=307)   470     21        0 (1447/1773-1487/1808)   0.0 W63-W73 N35-N37
      2 1402 1783 (N39W61 L=299)   449     63        0 (1386/1751-1430/1802)   6.3 W56-W65 N37-N39
      3 1717 1486 (N12W62 L=300)  1594    210        3 (1610/1447-1805/1534) -36.3 W53-W70 N07-N18
      4 1023 1355 (N26W09 L=247)   658    186        3 ( 931/1304-1092/1412) -52.1 W04-W11 N21-N31
      5 1063 1261 (N19W08 L=247)   419    180        0 ( 977/1232-1141/1281) -11.3 W02-W12 N17-N20
      6 1613 1147 (S02W38 L=276)   716    153        0 (1588/1066-1638/1207)  50.2 W35-W41 S05-N01
      7 1076 1072 (N08W04 L=242)   227    137        0 (1045/1015-1108/1136)  45.0 W01-W07 N05-N11
      8  323 1069 (N25E46 L=192)   365     46        1 ( 293/1042- 347/1096) -31.0 E49-E44 N22-N26
      9  264  909 (N16E54 L=184)   560    104        1 ( 226/ 862- 293/ 945) -59.0 E58-E52 N13-N18
     10 1861  871 (S26W56 L=294)   734    111        0 (1815/ 848-1895/ 907) -20.0 W49-W61 S28-S23
     11  467  847 (N09E38 L=200)   341     66        0 ( 442/ 810- 496/ 881) -68.2 E39-E36 N06-N11
     12  235  788 (N10E61 L=178)   409     36        0 ( 214/ 772- 264/ 806) -21.8 E63-E57 N08-N10
     13  438  537 (S08E51 L=187)   715    165        0 ( 384/ 473- 465/ 600) -90.0 E54-E49 S12-S03
     14 1383  391 (S42W06 L=244)  3577    731        4 (1085/ 221-1652/ 574)   7.4 E23-W31 S47-S36
     15 1190  520 (S27E03 L=235)   367     80        0 (1171/ 482-1210/ 566)  80.5 E04-E02 S29-S23
     16  638  472 (S16E39 L=199)  1164    236        3 ( 582/ 403- 707/ 545) -71.6 E43-E35 S21-S11
     17  782  428 (S22E32 L=207)   557    115        0 ( 731/ 397- 814/ 471)   7.1 E35-E27 S24-S19
     18  564  396 (S20E49 L=190)   464    150        1 ( 534/ 314- 591/ 446)  55.0 E54-E46 S25-S15
     19  762  237 (S36E46 L=192)   470     81        1 ( 749/ 194- 777/ 284)  38.7 E51-E41 S39-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