Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1195 1370 (N22W03 L=181)  1101    293        4 (1121/1241-1256/1472) -42.5 E04-W08 N15-N26
      2  973 1289 (N12E09 L=169)   243    137        0 ( 934/1229-1005/1344) -45.0 E12-E05 N08-N14
      3  199  998 (S21E58 L=120)   558    179        3 ( 167/ 961- 241/1044)  40.6 E61-E54 S23-S17
      4 1780  881 (N07W54 L=232)  1133    178        1 (1759/ 797-1802/ 961) -60.9 W51-W59 N02-N11
      5  544  871 (S21E27 L=151)   395    176        1 ( 490/ 831- 587/ 949)  60.9 E30-E24 S24-S15
      6 1313  790 (S09W22 L=200)   487    196        2 (1262/ 710-1344/ 855)  90.0 W18-W24 S15-S05
      7 1665  680 (S08W51 L=229)   600    139        2 (1610/ 652-1705/ 707)  12.5 W46-W55 S09-S05
      8 1121  573 (S27W17 L=195)   468    147        0 (1053/ 534-1168/ 635)  54.5 W13-W19 S30-S22
      9 1643  603 (S12W52 L=230)   661    174        4 (1590/ 568-1700/ 626)  24.0 W47-W56 S14-S09
     10 1500  587 (S17W41 L=219)   347     70        1 (1472/ 562-1526/ 620)  76.0 W39-W43 S18-S14
     11  596  597 (S38E20 L=158)   288    152        0 ( 534/ 572- 661/ 616)  26.6 E24-E14 S41-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