Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1311 1680 (N35W39 L=134)  1120    179        1 (1234/1655-1401/1707) -22.6 W32-W45 N31-N37
      2 1569 1603 (N23W54 L=149)   466     73        1 (1533/1575-1607/1630)   0.0 W49-W58 N22-N25
      3 1121 1343 (N19W13 L=107)  1369    401        2 ( 941/1281-1279/1451)   2.4 W01-W25 N16-N21
      4  584 1391 (N34E21 L= 73)  1700    301        2 ( 454/1309- 709/1446) -32.3 E32-E13 N25-N38
      5  155 1268 (N35E69 L= 26)   809     65        1 ( 134/1235- 186/1299) -20.0 E76-E62 N32-N36
      6 1745  888 (S21W45 L=140)  1760    367        1 (1674/ 759-1830/1043)  37.4 W37-W55 S28-S13
      7  875  723 (S12E15 L= 79)   645    194        2 ( 857/ 644- 900/ 827)  59.0 E17-E11 S15-S05
      8 1408  683 (S26W16 L=111)   271     22        0 (1392/ 658-1428/ 708) -90.0 W15-W17 S27-S23
      9 1342  528 (S34W09 L=104)   229      0       -1 (1317/ 510-1367/ 548) -26.6 W07-W11 S35-S32
     10  346  447 (S18E67 L= 27)   577     95        0 ( 318/ 419- 378/ 487) -90.0 E69-E64 S20-S14
  

    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