Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  677 1529 (N34E28 L= 11)   838    240        1 ( 583/1471- 750/1589) -23.6 E37-E21 N29-N37
      2 1233 1325 (N22W13 L= 52)  1486    363        3 (1062/1228-1399/1433) -26.6 W01-W23 N16-N28
      3  323 1347 (N19E53 L=346)   805    180        0 ( 303/1273- 342/1423) -63.4 E56-E49 N14-N24
      4  740 1391 (N25E21 L= 18)   357     56        1 ( 713/1369- 775/1416)  14.0 E22-E18 N23-N26
      5  911 1280 (N18E08 L= 31)   593    187        1 ( 857/1209- 986/1354) -45.0 E12-E03 N14-N22
      6  764 1239 (N15E17 L= 22)   489    157        1 ( 721/1161- 793/1294) -69.4 E20-E15 N10-N18
      7  517 1132 (N07E33 L=  6)   852    209        2 ( 483/1042- 576/1227) -55.0 E37-E28 N02-N12
      8  702 1066 (N04E20 L= 19)   200     73        1 ( 661/1050- 739/1086)  21.8 E22-E17 N02-N05
      9  427  969 (S03E40 L=360)   312     87        1 ( 390/ 953- 472/ 989)  -8.1 E42-E35 S04-S02
     10 1034  500 (S31W03 L= 42)   232     63        0 (1008/ 471-1053/ 525) -45.0 W00-W04 S33-S29
     11 1213  257 (S51W23 L= 62)   592    131        0 (1165/ 222-1247/ 285) -28.6 W17-W28 S55-S48
  

    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