Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1517 1776 (N35W69 L=355)   614    110        0 (1489/1740-1564/1790)  -9.5 W62-W77 N32-N35
      2  496 1663 (N58E32 L=255)  3303    453        5 ( 277/1577- 650/1749)   4.7 E75-E11 N55-N62
      3 1300 1549 (N27W32 L=318)   279     44        1 (1279/1531-1327/1571) -71.6 W29-W33 N25-N28
      4 1403 1470 (N20W35 L=321)   968    223        4 (1314/1420-1491/1518) -35.0 W29-W39 N14-N23
      5  791 1428 (N32E05 L=281)   271     65        0 ( 768/1399- 814/1458) -78.7 E06-E04 N29-N34
      6  365 1146 (N24E41 L=246)   254     93        0 ( 338/1123- 403/1186) -39.8 E43-E37 N21-N26
      7 1653 1013 (S12W39 L=325)   487    106        1 (1632/ 951-1675/1054)  59.0 W36-W40 S15-S09
      8  227  758 (N04E62 L=224)   607    100        1 ( 211/ 711- 250/ 798)  63.4 E64-E58 N00-N06
      9 1600  678 (S31W30 L=317)   282     76        0 (1558/ 662-1637/ 692)   0.0 W26-W33 S32-S30
     10  493  610 (S10E43 L=243)   328    104        1 ( 469/ 562- 519/ 655)  26.6 E47-E39 S12-S07
     11  995  400 (S34E19 L=268)  1355    392        0 ( 844/ 320-1124/ 533)  11.3 E30-E05 S37-S28
     12 1210  311 (S47E07 L=279)   358     43        1 (1180/ 292-1239/ 328) -11.3 E09-E04 S49-S45
  

    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