Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  861 1706 (N50W12 L=279)   425     53        1 ( 828/1684- 897/1722) -26.6 W09-W15 N49-N52
      2 1160 1424 (N23W19 L=286)   468     87        1 (1127/1380-1205/1462)  18.4 W16-W22 N20-N24
      3  954 1382 (N26W05 L=272)   850    316        3 ( 785/1353-1065/1408) -26.6 E05-W12 N22-N30
      4  245 1334 (N42E55 L=212)   830    108        1 ( 184/1311- 297/1353) -11.9 E66-E47 N39-N43
      5 1741  904 (S22W43 L=310)   858    158        0 (1706/ 830-1763/ 981)  58.0 W39-W45 S25-S17
      6 1507  899 (S15W25 L=292)   309    109        0 (1476/ 862-1549/ 944) -71.6 W23-W27 S17-S11
      7  341  784 (N06E49 L=218)   637     89        1 ( 307/ 746- 391/ 818)   0.0 E53-E45 N04-N08
      8  346  662 (S01E53 L=214)   675     75        2 ( 321/ 621- 374/ 695)  76.0 E55-E50 S03-N01
      9 1054  552 (S23E11 L=256)   599    195        0 ( 984/ 497-1098/ 637)  26.6 E16-E06 S25-S19
     10  927  518 (S22E19 L=248)   256     83        0 ( 893/ 496- 970/ 543) -53.1 E20-E17 S24-S19
     11 1284  438 (S38E00 L=267)   846    249        1 (1215/ 353-1352/ 527)  34.7 E07-W06 S42-S33
     12  745  392 (S26E36 L=231)   459    112        0 ( 730/ 329- 760/ 451)  56.3 E39-E33 S30-S21
  

    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