Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  488 1420 (N37E37 L=129)   292     56        0 ( 462/1403- 526/1433)  11.3 E39-E32 N35-N37
      2  873 1367 (N29E06 L=160)   260     46        0 ( 844/1348- 899/1384) -26.6 E07-E03 N27-N30
      3 1661 1205 (N07W45 L=211)   351     82        1 (1641/1171-1680/1252)  53.1 W43-W47 N05-N09
      4  534 1134 (N18E31 L=135)   240     74        0 ( 509/1094- 553/1163) -63.4 E32-E29 N15-N20
      5  451 1076 (N15E38 L=128)   543    185        0 ( 420/1002- 507/1154) -63.4 E40-E33 N10-N20
      6  313  919 (N06E51 L=115)   805    108        2 ( 277/ 878- 346/ 966)  35.5 E54-E47 N04-N09
      7 1602  858 (S13W36 L=201)   276     59        1 (1588/ 824-1619/ 888)  76.0 W34-W37 S14-S10
      8  316  520 (S20E65 L=100)   518     47        0 ( 301/ 494- 337/ 542)  18.4 E69-E60 S21-S17
      9 1591  394 (S48W43 L=208)  1356    146        2 (1526/ 362-1650/ 435)   8.1 W34-W50 S50-S44
     10  615  287 (S39E50 L=116)  1135    119        1 ( 560/ 255- 669/ 318)   9.5 E61-E42 S42-S36
  

    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