Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1593 1555 (N22W56 L= 36)  1526    258        4 (1548/1442-1633/1661)  31.4 W47-W65 N16-N27
      2 1390 1535 (N27W39 L= 18)  1171    187        1 (1314/1467-1454/1597) -74.7 W36-W41 N21-N32
      3 1036 1350 (N26W09 L=348)   381    148        0 ( 999/1281-1061/1400)  82.9 W07-W10 N21-N29
      4 1187 1357 (N22W18 L=358)   255     72        1 (1156/1330-1217/1388)  11.3 W16-W21 N20-N23
      5 1365 1326 (N16W28 L=  8)   628    245        0 (1271/1256-1455/1382) -57.5 W23-W32 N10-N21
      6  626 1094 (N20E23 L=316)   501    129        2 ( 601/1028- 657/1154)  66.8 E25-E20 N15-N22
      7  833  712 (S07E18 L=321)   266    129        0 ( 806/ 665- 885/ 757)  33.7 E20-E14 S09-S05
      8 1350  700 (S20W12 L=351)   552    248        1 (1224/ 674-1426/ 752)  -8.7 W04-W17 S21-S16
      9 1533  687 (S26W25 L=  4)   827    183        1 (1460/ 633-1614/ 725)   4.4 W18-W31 S27-S23
     10  745  281 (S34E43 L=296)  2273    384        1 ( 709/ 143- 797/ 476)  33.1 E71-E28 S48-S20
  

    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