Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  742 1456 (N39E10 L=119)   580    109        2 ( 695/1420- 786/1483) -29.7 E13-E06 N36-N41
      2  489 1236 (N30E32 L= 97)  1056    255        1 ( 421/1153- 560/1330) -54.5 E37-E26 N23-N37
      3 1025 1248 (N20W05 L=134)   547    200        1 ( 937/1202-1126/1294) -38.7 W00-W10 N15-N24
      4 1309 1054 (N03W17 L=146)   982    234        1 (1203/ 963-1376/1163) -70.3 W13-W20 S03-N11
      5  508 1035 (N17E32 L= 96)   304    121        0 ( 482/ 997- 546/1094) -74.1 E33-E30 N14-N21
      6 1848 1040 (S13W58 L=186)   796     79        2 (1820/1005-1867/1080)  76.0 W54-W60 S14-S10
      7 1307  721 (S17W11 L=139)   753    158        1 (1236/ 654-1344/ 772) -54.5 W06-W12 S21-S14
      8  305  514 (S12E66 L= 63)  1267    260        1 ( 282/ 462- 343/ 576)  33.7 E73-E60 S16-S07
      9 1048  359 (S36E15 L=113)   854    165        2 (1013/ 304-1082/ 424)  45.0 E20-E11 S39-S30
     10  604  328 (S29E48 L= 81)   945    141        2 ( 543/ 294- 649/ 358) -29.7 E52-E43 S31-S25
  

    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