Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  578 1735 (N42E44 L=184)   445     51        2 ( 550/1711- 601/1757)  45.0 E46-E40 N40-N44
      2 1424 1645 (N38W29 L=256)  2282    392        4 (1276/1537-1589/1708) -22.6 W16-W41 N31-N43
      3  330 1508 (N25E56 L=172)   881    113        0 ( 310/1457- 352/1566) -45.0 E61-E52 N21-N29
      4  900 1536 (N27E11 L=217)   266     68        1 ( 867/1515- 938/1559) -21.8 E12-E07 N26-N29
      5 1682 1418 (N24W46 L=274)   607    116        0 (1661/1371-1714/1477) -71.6 W43-W48 N21-N28
      6 1519 1171 (N07W31 L=258)   454     84        0 (1488/1134-1538/1219)  80.5 W28-W31 N04-N10
      7  405 1092 (S01E40 L=187)   244     70        1 ( 385/1052- 420/1130) -71.6 E41-E39 S03-N00
      8  255  877 (S14E55 L=173)  1433    239        4 ( 213/ 788- 301/ 981)  67.4 E61-E50 S19-S07
      9  491  923 (S12E34 L=194)   394    106        0 ( 458/ 871- 525/ 959)  45.0 E36-E31 S14-S09
     10 1618  847 (S12W40 L=268)   317     44        0 (1585/ 827-1651/ 870)   0.0 W37-W42 S12-S10
     11  910  782 (S19E06 L=222)   348    112        1 ( 871/ 753- 962/ 825)  38.7 E08-E03 S20-S16
     12 1717  763 (S16W50 L=278)   520    101        3 (1695/ 721-1745/ 803) -45.0 W47-W53 S18-S13
     13  454  630 (S30E41 L=187)  1009    136        2 ( 419/ 563- 493/ 706)  84.3 E42-E37 S34-S24
     14 1716  633 (S24W54 L=282)   353     48        1 (1696/ 605-1734/ 657) -71.6 W52-W56 S25-S22
     15  208  606 (S31E73 L=155)  1217     89        3 ( 185/ 571- 235/ 644)  11.3 E85-E65 S32-S28
     16 1364  593 (S29W26 L=253)   600    112        2 (1320/ 558-1396/ 640)  63.4 W22-W28 S31-S25
     17 1111  588 (S31W08 L=236)   560    187        0 (1017/ 572-1198/ 606)   4.8 W01-W13 S31-S29
     18  858  496 (S38E10 L=218)   280     67        0 ( 825/ 467- 883/ 529)  51.3 E12-E08 S40-S35
     19  572  279 (S56E50 L=178)   469     39        1 ( 548/ 261- 603/ 303)   8.7 E56-E43 S58-S54
  

    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