Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1451 1572 (N37W38 L= 31)   428     73        2 (1427/1537-1477/1604)  45.0 W36-W41 N34-N38
      2 1385 1486 (N31W29 L= 22)   607    130        1 (1329/1456-1423/1526)  -7.1 W24-W32 N29-N34
      3 1439 1367 (N23W30 L= 23)   259    113        0 (1413/1329-1473/1417) -71.6 W29-W32 N20-N26
      4 1602 1321 (N19W42 L= 35)   398    116        1 (1569/1290-1653/1356) -26.6 W40-W46 N16-N21
      5  288 1096 (N11E52 L=300)   552    179        0 ( 257/1023- 349/1143)  41.2 E54-E46 N06-N13
      6  388 1008 (N05E43 L=310)  1045    212        2 ( 365/ 916- 416/1113) -84.8 E44-E40 N00-N11
      7  114 1014 (N05E76 L=276)  2980    178        1 (  91/ 924- 142/1087) -48.8 E84-E70 N00-N09
      8  863  950 (N01E10 L=342)   451    218        0 ( 806/ 870- 885/1039)  60.9 E14-E08 S03-N06
      9 1906  854 (S13W71 L= 64)   723     59        2 (1881/ 830-1932/ 877)  11.3 W66-W77 S14-S12
     10  260  749 (S11E58 L=294)  2112    346        4 ( 131/ 640- 378/ 818) -17.2 E77-E48 S17-S06
     11 1049  699 (S16E00 L=352)   197     63        1 (1035/ 660-1067/ 723)  76.0 E01-W01 S18-S14
     12  319  597 (S21E57 L=295)   970    151        2 ( 254/ 558- 390/ 637) -17.1 E63-E50 S23-S18
  

    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