Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  897 1773 (N37E35 L= 16)  3037    583        3 ( 851/1652- 969/1930) -38.2 E54-E24 N28-N50
      2 1599 1514 (N42W25 L= 76)   303     63        1 (1575/1485-1623/1534) -45.0 W22-W27 N40-N44
      3 1772 1388 (N40W43 L= 94)  2228    400        5 (1645/1310-1859/1474)  -4.8 W29-W53 N36-N42
      4 1344 1414 (N26W08 L= 59)   239    112        0 (1316/1350-1363/1456) -56.3 W05-W09 N21-N28
      5  773 1351 (N06E23 L= 28)   411     67        2 ( 747/1315- 798/1385)  76.0 E24-E21 N04-N08
      6 1844 1258 (N34W51 L=102)   749    136        0 (1801/1227-1870/1285)  32.0 W46-W54 N31-N36
      7 1042 1158 (N03E03 L= 49)   195     95        0 ( 999/1147-1097/1169)  21.8 E05-W00 N01-N04
      8 1121 1147 (N04W02 L= 53)   209     12        0 (1102/1126-1139/1169) -71.6 W01-W02 N02-N05
      9 1571 1038 (N11W32 L= 83)   502    124        1 (1530/ 987-1602/1089) -38.7 W27-W34 N07-N13
     10 1660 1040 (N14W38 L= 89)   878    237        1 (1623/ 965-1714/1089)  82.9 W35-W42 N09-N17
     11  827  952 (S14E09 L= 42)   279     73        2 ( 810/ 906- 850/ 990) -78.7 E10-E08 S17-S12
     12 1690  904 (N07W44 L= 95)   881    144        1 (1642/ 848-1732/ 964)  76.0 W41-W48 N03-N11
     13 1829  904 (N12W57 L=108)   602     99        0 (1804/ 862-1856/ 946)  90.0 W54-W58 N09-N15
     14  676  858 (S24E16 L= 35)   696    218        0 ( 578/ 825- 776/ 906)  36.9 E22-E10 S27-S18
     15  344  875 (S31E41 L= 10)   334    111        1 ( 308/ 848- 384/ 901)  33.7 E44-E38 S33-S28
     16  853  705 (S28E01 L= 50)   452     68        2 ( 813/ 673- 884/ 737)  68.2 E03-W01 S30-S25
     17  636  631 (S39E14 L= 37)   439     61        0 ( 606/ 611- 683/ 658)  -8.1 E16-E09 S40-S37
     18 1176  302 (S42W39 L= 90)   297     59        0 (1150/ 279-1199/ 320)   0.0 W36-W41 S43-S41
  

    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