Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  675 1616 (N50E07 L=265)   346     66        2 ( 658/1578- 693/1650)  59.0 E08-E05 N47-N53
      2 1356 1610 (N28W40 L=312)   467     59        2 (1323/1583-1393/1631) -45.0 W38-W43 N26-N30
      3  446 1163 (N26E33 L=238)   252      6        0 ( 427/1144- 465/1181) -45.0 E35-E31 N25-N27
      4  518 1088 (N20E29 L=243)   448     95        0 ( 473/1063- 575/1105)  -7.1 E32-E24 N19-N21
      5  144  969 (N21E68 L=204)   567     81        1 ( 109/ 936- 169/1007) -29.1 E74-E64 N18-N23
      6  219  930 (N17E58 L=214)   403     54        1 ( 199/ 902- 241/ 957)  36.9 E60-E55 N15-N18
      7 1876  796 (S34W60 L=331)  1371    126        2 (1849/ 735-1901/ 852)  66.8 W55-W63 S37-S30
      8 1713  646 (S38W39 L=311)   349     32        0 (1696/ 620-1728/ 672)  76.0 W36-W40 S39-S35
      9 1217  549 (S28E01 L=271)   418    164        0 (1164/ 492-1281/ 599)  26.6 E05-W04 S30-S25
     10 1352  234 (S59E03 L=269)   538     53        1 (1320/ 218-1382/ 255) -50.2 E06-W00 S62-S56
  

    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