Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1376 1631 (N42W35 L=289)   299     44        0 (1346/1618-1403/1647) -11.3 W31-W37 N41-N43
      2 1772 1272 (N15W58 L=312)  1471    282        2 (1674/1179-1841/1359) -38.2 W50-W64 N08-N20
      3  323 1161 (N14E50 L=205)   307     71        0 ( 287/1146- 352/1174)   0.0 E53-E47 N13-N15
      4 1100  865 (S06W04 L=259)   295    120        0 (1038/ 832-1152/ 887) -23.2 W00-W07 S07-S04
      5 1445  640 (S22W27 L=282)   486     65        1 (1418/ 604-1478/ 669)  63.4 W25-W29 S24-S20
      6  582  593 (S22E33 L=222)   380     89        1 ( 551/ 558- 615/ 623) -59.0 E34-E30 S24-S19
      7 1308  524 (S29W18 L=272)   568    134        3 (1258/ 482-1363/ 574)  39.8 W14-W21 S32-S26
      8  955  499 (S29E07 L=247)   678    278        0 ( 852/ 437-1072/ 556)  26.6 E15-W01 S33-S25
      9  559  517 (S27E37 L=218)   340    110        0 ( 535/ 469- 598/ 554)  40.6 E39-E32 S30-S24
     10 1580  393 (S43W49 L=303)   591     46        0 (1556/ 372-1608/ 417) -20.0 W43-W54 S45-S40
  

    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