Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1151 1756 (N54W01 L= 83)   584    109        1 (1096/1734-1207/1775) -10.3 E05-W06 N52-N55
      2  434 1493 (N24E51 L= 31)  1994    285        4 ( 355/1397- 479/1626) -51.3 E57-E44 N18-N33
      3 1281 1430 (N30W13 L= 95)  1556    468        1 (1098/1282-1397/1551) -35.5 W00-W22 N20-N36
      4 1766 1392 (N32W58 L=140)   444     45        0 (1747/1364-1783/1420) -71.6 W56-W60 N29-N33
      5 1102 1236 (N15W03 L= 85)   368    149        0 (1071/1190-1149/1301) -50.2 W00-W05 N12-N18
      6 1543 1210 (N18W33 L=115)   334    170        0 (1494/1180-1570/1255) -38.7 W29-W35 N15-N20
      7 1219  683 (S18W16 L= 98)   739    178        0 (1145/ 626-1293/ 738)  45.0 W12-W20 S22-S13
      8 1697  655 (S15W53 L=135)   473     70        2 (1660/ 627-1726/ 681)   0.0 W49-W56 S17-S13
      9 1595  579 (S21W47 L=129)   403     39        0 (1575/ 551-1609/ 611)  90.0 W45-W48 S23-S18
     10 1659  547 (S23W55 L=137)   376     42        1 (1641/ 520-1680/ 575) -63.4 W52-W57 S24-S20
     11  404  513 (S39E48 L= 34)  1029    166        2 ( 348/ 482- 474/ 556)  20.2 E57-E37 S42-S35
  

    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