Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1435 1559 (N41W23 L=157)   583    151        0 (1386/1518-1484/1597)  33.7 W18-W27 N37-N44
      2  413 1448 (N17E50 L= 84)   543    125        0 ( 370/1418- 448/1499)  63.4 E53-E47 N14-N20
      3 1338 1460 (N32W15 L=149)   286     43        1 (1316/1431-1358/1487) -71.6 W13-W16 N30-N34
      4  775 1341 (N16E20 L=114)   225     93        0 ( 724/1326- 815/1357)  18.4 E23-E17 N13-N17
      5  992 1248 (N13E05 L=129)   338     75        0 ( 944/1227-1033/1273)  21.8 E07-E02 N11-N14
      6 1000 1150 (N07E03 L=131)   779    294        2 ( 890/1038-1099/1229) -31.0 E11-W04 N01-N10
      7  641  940 (S11E23 L=112)   347    129        0 ( 577/ 922- 695/ 963)  15.9 E26-E19 S11-S08
      8 1001  852 (S11W01 L=136)   725    170        4 ( 972/ 767-1029/ 930) -78.7 E01-W04 S15-S05
      9  452  566 (S38E37 L= 98)   458    110        0 ( 412/ 539- 491/ 597)  32.0 E40-E32 S40-S35
  

    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