Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  235 1318 (N41E62 L= 42)   427     53        3 ( 205/1300- 264/1330) -10.3 E67-E56 N39-N41
      2  366 1285 (N36E43 L= 60)   327     90        0 ( 331/1259- 410/1301)   7.1 E47-E38 N35-N37
      3  383 1098 (N24E42 L= 62)   402     98        1 ( 352/1055- 407/1153) -81.9 E44-E39 N21-N28
      4 1707  989 (S12W43 L=146)   334     49        1 (1677/ 965-1738/1004) -21.8 W40-W45 S14-S11
      5  164  896 (N15E69 L= 34)  1126    164        3 ( 140/ 821- 199/ 972) -83.7 E72-E65 N10-N19
      6  417  839 (N08E43 L= 60)   497     94        1 ( 391/ 807- 452/ 874)  18.4 E45-E39 N06-N10
      7 1198  620 (S21W01 L=104)  1618    542        3 (1026/ 471-1338/ 786)  27.5 E13-W12 S26-S13
      8  901  445 (S25E23 L= 81)  1540    402        2 ( 838/ 373-1018/ 604)  56.3 E27-E15 S31-S15
      9 1248  433 (S35E01 L=103)   300     79        1 (1216/ 407-1278/ 460)  21.8 E03-W02 S36-S33
     10  587  265 (S32E60 L= 43)   924     77        2 ( 565/ 225- 610/ 305)  45.0 E65-E55 S36-S28
  

    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