Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1363 1760 (N43W50 L= 82)   515    117        1 (1342/1724-1392/1782)  15.9 W46-W53 N41-N44
      2 1270 1418 (N25W25 L= 56)   404     86        2 (1233/1385-1305/1447)  18.4 W21-W27 N23-N26
      3  899 1388 (N31E00 L= 31)   305    101        0 ( 850/1373- 959/1398) -26.6 E03-W03 N29-N32
      4  371 1355 (N39E43 L=349)   987    197        2 ( 310/1310- 443/1387) -17.1 E48-E35 N34-N41
      5  402 1108 (N22E40 L=352)   574    216        1 ( 311/1092- 478/1129)  -8.1 E47-E33 N20-N24
      6  716  994 (N11E19 L= 13)   375     71        1 ( 671/ 977- 754/1015) -31.0 E21-E16 N09-N12
      7 1461  391 (S45W19 L= 51)   743    186        1 (1389/ 375-1530/ 411) -29.7 W12-W26 S49-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