Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  603 1540 (N34E35 L=331)   507    126        1 ( 587/1479- 624/1601) -48.4 E39-E31 N29-N38
      2 1593 1480 (N33W44 L= 50)   863    139        2 (1524/1439-1643/1519) -24.0 W38-W48 N30-N36
      3  886 1480 (N31E12 L=354)   492    201        1 ( 803/1421- 960/1511) -18.4 E18-E06 N27-N33
      4 1324 1433 (N30W20 L= 26)   399    195        0 (1229/1387-1399/1462) -18.4 W12-W24 N26-N31
      5  460 1310 (N18E41 L=326)   302     97        0 ( 434/1273- 496/1348) -45.0 E43-E36 N15-N20
      6 1119 1199 (N14W05 L= 12)   308     75        1 (1088/1179-1160/1229) -26.6 W03-W07 N12-N15
      7 1793 1040 (N05W55 L= 62)   421     95        0 (1752/1001-1830/1068)  20.6 W51-W59 N02-N07
      8  470  896 (S08E36 L=330)   825    213        2 ( 412/ 799- 510/ 966)  65.6 E40-E32 S14-S03
      9 1692  906 (S03W46 L= 52)   820    111        3 (1639/ 867-1743/ 948) -16.7 W41-W51 S05-N00
     10 1064  800 (S11W03 L= 10)   257    104        0 (1026/ 763-1095/ 827) -38.7 W00-W05 S13-S09
     11 1819  768 (S12W62 L= 68)   490     63        2 (1800/ 741-1845/ 798) -53.1 W59-W66 S13-S09
     12  873  672 (S20E09 L=358)   890    218        5 ( 812/ 593- 922/ 765) -76.0 E12-E05 S25-S13
     13  981  693 (S18E02 L=  5)   622    116        3 ( 940/ 642-1028/ 734) -50.2 E04-W01 S21-S15
     14  424  564 (S30E45 L=321)  1609    276        2 ( 368/ 507- 539/ 657) -28.3 E49-E34 S34-S23
     15 1166  485 (S32W13 L= 19)  4253    700        8 ( 941/ 382-1393/ 647)  28.1 E03-W27 S40-S19
     16  759  578 (S27E17 L=350)   741    165        1 ( 699/ 530- 831/ 617)  35.0 E21-E11 S30-S23
     17 1512  496 (S30W40 L= 46)  1503    307        3 (1390/ 431-1639/ 555)  26.6 W30-W49 S35-S25
  

    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