Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  647 1613 (N20E42 L=179)   368     78        2 ( 615/1583- 678/1634)  -8.1 E45-E38 N18-N21
      2 1320 1605 (N39W02 L=223)   475    100        0 (1265/1589-1375/1619)  20.6 E02-W06 N36-N41
      3  901 1561 (N23E23 L=198)   270     90        0 ( 889/1512- 913/1602) -45.0 E25-E20 N20-N25
      4 1251 1434 (N25W02 L=223)   366    106        0 (1219/1388-1287/1475)  81.9 W01-W03 N21-N28
      5  669 1131 (S07E23 L=198)   207     61        1 ( 642/1101- 698/1158) -26.6 E25-E20 S08-S06
      6  214 1062 (S21E59 L=162)  1158    172        2 ( 196/ 979- 243/1144) -78.7 E61-E56 S26-S16
      7  440  903 (S26E34 L=186)   462    127        0 ( 391/ 869- 497/ 939)  -5.7 E39-E29 S27-S24
      8  738  749 (S28E10 L=211)   359     94        1 ( 689/ 724- 784/ 783)   0.0 E13-E06 S28-S26
      9 1260  719 (S15W22 L=243)   242     47        3 (1238/ 691-1282/ 743) -71.6 W20-W23 S16-S13
     10  307  579 (S52E52 L=169)   561     52        2 ( 290/ 551- 327/ 616)  68.2 E56-E48 S53-S48
     11 1400  457 (S26W43 L=264)   325     72        0 (1372/ 424-1432/ 495) -80.5 W42-W44 S28-S22
  

    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