Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  823 1812 (N38E45 L=217)  1383    181        2 ( 750/1784- 913/1845)   3.2 E54-E35 N36-N40
      2  476 1693 (N21E63 L=198)   520     67        1 ( 459/1671- 498/1723) -20.6 E67-E58 N19-N23
      3  729 1656 (N24E39 L=223)   289     58        1 ( 705/1632- 761/1682)  76.0 E40-E37 N22-N26
      4  290 1542 (N08E72 L=189)  1545    150        3 ( 244/1473- 331/1584)  69.4 E77-E67 N02-N11
      5 1270 1376 (N22W05 L=266)   497    179        0 (1232/1335-1332/1433) -23.2 W01-W09 N20-N25
      6  193 1300 (S08E69 L=192)   742     63        2 ( 164/1281- 226/1327)  20.6 E74-E64 S09-S06
      7  243  969 (S27E54 L=207)   623     86        1 ( 218/ 945- 279/1013) -26.6 E57-E50 S28-S24
      8 1804  840 (N08W57 L=319)  1230    210        4 (1749/ 770-1852/ 901)  60.9 W53-W63 N02-N11
      9  711  577 (S40E07 L=254)   371    107        0 ( 686/ 531- 754/ 624) -29.7 E09-E02 S41-S36
     10 1342  540 (S23W35 L=296)   351    107        0 (1298/ 521-1393/ 568)   7.1 W30-W38 S23-S20
  

    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