Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1236 1748 (N46W07 L=260)  1212    167        2 (1155/1710-1311/1777)   0.0 W00-W13 N43-N48
      2  353 1535 (N21E57 L=196)  1266    125        2 ( 329/1473- 390/1605) -45.0 E64-E52 N16-N25
      3 1012 1437 (N19E05 L=247)   955    286        1 ( 908/1335-1119/1528) -33.7 E13-W02 N13-N24
      4 1467 1427 (N25W25 L=278)   394    103        1 (1429/1389-1508/1471) -18.4 W21-W27 N22-N27
      5 1723 1331 (N23W46 L=299)  1445    303        5 (1687/1188-1763/1464)  86.8 W42-W50 N13-N31
      6 1654 1392 (N26W40 L=293)   349     59        0 (1630/1366-1683/1422) -45.0 W37-W42 N24-N27
      7 1809 1332 (N25W55 L=308)   509     70        0 (1790/1293-1836/1357)  39.8 W52-W59 N21-N26
      8 1848 1141 (N13W58 L=311)  1147    226        1 (1801/1045-1902/1242) -42.3 W52-W65 N08-N18
      9  329  679 (S33E51 L=202)  1264    313        0 ( 275/ 532- 391/ 804)  36.9 E67-E42 S42-S24
     10 1826  660 (S15W68 L=321)  3490    198        6 (1760/ 572-1882/ 748) -84.8 W57-W79 S20-S09
     11 1176  614 (S28W15 L=268)  1279    401        1 (1036/ 502-1336/ 721)  48.4 W07-W23 S37-S19
     12  854  463 (S43E05 L=247)   227     43        0 ( 842/ 434- 871/ 491)  90.0 E06-E04 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