Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  530 1700 (N27E58 L=100)   755    109        2 ( 504/1650- 551/1746) -24.8 E65-E51 N24-N30
      2  345 1587 (N16E69 L= 89)  2693    123        4 ( 302/1519- 390/1637) -56.3 E78-E63 N11-N20
      3  479 1564 (N17E52 L=106)  1277    157        2 ( 423/1519- 546/1597)  49.4 E57-E47 N13-N20
      4  556 1518 (N16E43 L=115)   453     64        2 ( 537/1486- 580/1549)  76.0 E45-E40 N14-N18
      5  483 1487 (N13E47 L=111)   371     62        0 ( 459/1459- 502/1519)  78.7 E48-E45 N10-N15
      6 1675 1460 (N41W40 L=198)   350     29        1 (1654/1445-1696/1479)  45.0 W36-W42 N39-N43
      7 1789 1296 (N33W51 L=210)   521     52        0 (1762/1272-1811/1316) -14.0 W47-W54 N31-N34
      8  242 1217 (S08E60 L= 99)  1101    130        1 ( 197/1177- 297/1251)  10.3 E64-E53 S10-S05
      9  469 1014 (S15E35 L=124)   542    188        0 ( 417/ 952- 517/1073)  60.9 E37-E32 S19-S10
     10  222  916 (S26E57 L=101)   775     91        1 ( 193/ 872- 257/ 959)  45.0 E61-E53 S29-S23
     11 1324  730 (S12W25 L=183)   284     96        1 (1302/ 691-1360/ 762)  71.6 W23-W26 S15-S09
     12  998  671 (S24W07 L=165)   199     14        0 ( 976/ 656-1019/ 688)   0.0 W05-W08 S24-S22
     13  903  621 (S29W02 L=161)   262     79        1 ( 866/ 592- 938/ 650) -21.8 W00-W05 S30-S28
     14  412  559 (S47E39 L=120)  1371    244        3 ( 336/ 527- 531/ 603)  16.3 E49-E25 S49-S41
     15 1427  559 (S20W38 L=197)   502     72        1 (1401/ 522-1451/ 596)  90.0 W36-W39 S22-S17
     16  909  499 (S37W06 L=165)   252     36        1 ( 886/ 479- 931/ 522)  56.3 W05-W07 S38-S35
     17  879  284 (S55W14 L=173)  1705    232        2 ( 834/ 198- 923/ 397) -64.5 W10-W21 S65-S44
  

    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