Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  411 1581 (N17E60 L=109)  4306    425       10 ( 291/1486- 514/1706)  86.6 E70-E50 N09-N26
      2  508 1650 (N23E56 L=114)   808     96        1 ( 468/1615- 534/1685)  90.0 E59-E52 N20-N26
      3 1637 1652 (N58W43 L=213)   699     61        1 (1615/1612-1661/1676)  32.7 W38-W52 N53-N62
      4  570 1491 (N15E41 L=128)   354     65        2 ( 547/1461- 599/1511)  90.0 E43-E39 N12-N16
      5 1705 1312 (N32W41 L=211)   869    186        2 (1665/1245-1762/1377) -35.0 W37-W47 N28-N35
      6  171 1237 (S08E70 L= 99)   795     89        0 ( 129/1214- 204/1258)   0.0 E79-E64 S09-S06
      7  345 1092 (S13E46 L=123)   285     77        0 ( 330/1052- 365/1122)  63.4 E47-E44 S15-S11
      8  144  935 (S27E70 L= 99)   682     62        1 ( 126/ 899- 170/ 967)  39.8 E75-E65 S29-S24
      9  854  721 (S24E03 L=166)   360    134        0 ( 831/ 676- 880/ 769) -59.0 E04-E01 S26-S21
     10  768  674 (S29E08 L=162)   227     88        0 ( 726/ 653- 811/ 695)   0.0 E10-E04 S29-S28
     11 1427  633 (S15W35 L=205)   296     56        0 (1397/ 608-1454/ 657)  63.4 W33-W36 S17-S13
     12 1294  625 (S19W27 L=196)   246     43        0 (1280/ 597-1308/ 653) -71.6 W25-W27 S20-S17
     13  314  589 (S47E52 L=118)  1323    142        1 ( 261/ 562- 378/ 631)   0.0 E61-E42 S49-S44
     14  744  542 (S39E07 L=163)   282     56        0 ( 727/ 514- 767/ 566)  78.7 E07-E04 S41-S36
     15  765  384 (S50E01 L=169)   725    105        2 ( 740/ 335- 805/ 424) -54.5 E03-W02 S54-S46
     16 1147  401 (S37W27 L=197)   406     66        1 (1104/ 385-1187/ 417)   8.1 W23-W30 S38-S35
     17 1537  368 (S28W61 L=230)   982    133        1 (1491/ 331-1595/ 415)  82.9 W56-W63 S31-S23
     18  789  284 (S59W06 L=176)   786    124        0 ( 737/ 250- 847/ 314)  31.0 W00-W10 S62-S53
  

    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