Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  962 1539 (N26E08 L=197)   637    258        0 ( 863/1461-1053/1590) -24.8 E14-E01 N21-N29
      2  295 1375 (N13E54 L=151)  2961    367       13 ( 219/1242- 391/1526) -42.0 E66-E42 N04-N23
      3  315 1494 (N21E57 L=148)   383     36        1 ( 298/1470- 329/1520) -71.6 E58-E54 N19-N22
      4  553 1444 (N18E34 L=171)   618    180        2 ( 492/1378- 599/1515)  60.9 E37-E31 N13-N22
      5  470 1417 (N16E40 L=166)   419     71        2 ( 440/1385- 502/1443)  63.4 E42-E37 N13-N17
      6  390 1399 (N15E46 L=159)   435     99        3 ( 363/1351- 413/1441) -80.5 E47-E43 N11-N17
      7 1516 1283 (N15W29 L=234)   332     53        0 (1484/1254-1539/1307) -14.0 W26-W30 N12-N15
      8  221 1076 (S05E56 L=149)   419    115        1 ( 181/1028- 265/1109)  26.6 E60-E51 S07-S03
      9  628 1046 (S06E24 L=181)   343    119        2 ( 591/ 996- 654/1087) -59.0 E26-E22 S08-S03
     10  195  906 (S15E60 L=145)  1643    256        2 ( 155/ 808- 247/1030)  52.4 E66-E53 S20-S07
     11  548  974 (S11E29 L=176)   233     57        1 ( 529/ 946- 575/ 994)  90.0 E30-E27 S12-S09
     12  644  933 (S12E23 L=182)   348     94        1 ( 604/ 899- 684/ 969)  38.7 E25-E20 S14-S10
     13  406  895 (S16E40 L=165)  1166    326        0 ( 331/ 805- 471/ 948)  35.0 E46-E34 S21-S12
     14  255  877 (S17E54 L=152)   643    111        2 ( 218/ 844- 284/ 933)  39.8 E57-E50 S18-S13
     15  552  829 (S19E29 L=176)   310    107        1 ( 512/ 800- 606/ 866) -26.6 E31-E25 S20-S16
     16  216  674 (S29E64 L=141)   548     38        1 ( 197/ 645- 232/ 699)  38.7 E67-E61 S31-S27
     17 1356  504 (S34W29 L=234)   526    161        2 (1281/ 466-1415/ 538) -14.0 W22-W34 S35-S32
     18 1389  238 (S53W51 L=256)   791    104        3 (1338/ 218-1444/ 256)   3.6 W42-W58 S54-S51
  

    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