Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1504 1771 (N35W58 L=175)   982    220        1 (1474/1728-1535/1810)  15.3 W52-W63 N33-N37
      2  404 1579 (N52E39 L= 77)   404     58        0 ( 376/1560- 429/1596) -20.6 E43-E35 N49-N53
      3  418 1406 (N38E35 L= 82)   289     44        0 ( 396/1379- 445/1426)  26.6 E36-E32 N36-N39
      4 1014 1330 (N20W06 L=123)   354    162        0 ( 936/1304-1094/1351) -33.7 W01-W10 N16-N22
      5 1830 1224 (S04W60 L=176)  3141    642        7 (1726/1024-1904/1346)  20.6 W44-W76 S12-N01
      6 1563 1290 (N06W39 L=156)   557    114        1 (1517/1233-1604/1339) -74.1 W36-W40 N02-N09
      7  897 1203 (N15E04 L=113)   184     81        1 ( 853/1189- 944/1221) -31.0 E06-E01 N13-N16
      8  180 1015 (N18E59 L= 57)  2497    539        8 (  97/ 919- 307/1142) -41.6 E68-E48 N10-N26
      9  960  504 (S27E16 L=101)   775    217        2 ( 919/ 428- 992/ 616)  54.0 E20-E11 S31-S20
     10 1620  534 (S42W33 L=149)   316     45        1 (1596/ 514-1642/ 557)  14.0 W30-W34 S43-S40
     11  467  358 (S26E58 L= 59)  1556    287        6 ( 434/ 285- 493/ 422)  39.3 E63-E51 S31-S21
     12  859  233 (S44E37 L= 80)   388     81        0 ( 845/ 206- 874/ 267)  51.3 E39-E34 S46-S41
  

    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