Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  401 1510 (N31E54 L= 37)   503     84        2 ( 364/1489- 444/1526)  -5.7 E58-E48 N29-N31
      2  743 1374 (N24E20 L= 70)   529    174        1 ( 661/1335- 820/1407)  -4.8 E26-E14 N21-N25
      3  884 1249 (N16E10 L= 81)   233    117        1 ( 822/1226- 932/1265)  -8.1 E13-E06 N15-N17
      4 1325 1198 (N14W19 L=109)   250     96        0 (1284/1178-1364/1216) -11.3 W16-W21 N13-N15
      5 1232  844 (S08W13 L=104)   199     93        0 (1188/ 821-1264/ 873)  31.0 W10-W15 S09-S05
      6 1825  728 (S15W65 L=155)   670    123        1 (1797/ 687-1856/ 766)  51.3 W61-W68 S17-S12
      7  928  466 (S34E05 L= 85)   974    296        1 ( 880/ 353- 998/ 582) -71.6 E08-W00 S43-S25
      8 1519  533 (S28W39 L=129)   385    100        2 (1482/ 499-1550/ 580)  63.4 W36-W40 S30-S24
  

    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