Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1284 1780 (N42W52 L=170)   369     58        0 (1271/1750-1300/1809)  20.6 W48-W57 N39-N43
      2  174 1036 (N26E65 L= 54)  4841    660        9 ( 117/ 803- 227/1319) -84.6 E75-E59 N12-N44
      3  357 1042 (N23E43 L= 75)  1230    304        3 ( 248/ 922- 444/1118) -49.1 E54-E37 N14-N29
      4  240  950 (N20E57 L= 62)  1278    232        3 ( 195/ 845- 279/1032) -76.0 E61-E53 N13-N25
      5 1497  900 (S13W24 L=142)  1105    319        0 (1373/ 837-1617/1003)   5.7 W14-W34 S15-S10
      6 1747  921 (S20W43 L=162)   842    172        0 (1704/ 849-1792/1000)  33.7 W38-W49 S23-S17
      7 1895  933 (S26W62 L=180)   706     63        0 (1872/ 908-1929/ 961) -29.1 W58-W68 S28-S23
      8 1794  700 (S37W48 L=166)   440     51        1 (1767/ 678-1820/ 722)  11.3 W44-W51 S38-S35
      9  858  390 (S26E30 L= 89)  3434    897        0 ( 684/ 250-1138/ 621)  11.8 E47-E04 S34-S18
     10 1666  565 (S41W33 L=151)   693    123        0 (1617/ 537-1721/ 606)   4.8 W27-W39 S42-S39
     11 1193  475 (S30E06 L=113)   453    155        0 (1131/ 436-1234/ 533)  24.0 E10-E01 S31-S27
  

    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