Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  584 1732 (N41E43 L=296)   403     44        0 ( 556/1711- 612/1750) -20.6 E47-E39 N39-N42
      2  483 1707 (N39E53 L=286)   571     59        0 ( 463/1668- 508/1734) -80.5 E55-E50 N35-N41
      3 1550 1566 (N34W37 L= 17)   350     50        0 (1534/1534-1565/1596) -90.0 W35-W38 N31-N36
      4  326 1549 (N27E60 L=280)   599     47        2 ( 295/1525- 348/1568) -16.7 E65-E55 N25-N28
      5  314  842 (S17E49 L=290)   449    101        0 ( 302/ 789- 328/ 893)  63.4 E50-E47 S20-S14
      6  208  761 (S22E63 L=277)   921     77        1 ( 161/ 739- 256/ 787)   0.0 E71-E56 S23-S20
      7 1715  644 (S22W54 L= 34)   325     22        0 (1696/ 619-1731/ 664) -33.7 W51-W56 S23-S21
      8 1665  526 (S30W55 L= 34)  1087     74        3 (1631/ 494-1705/ 558) -10.3 W50-W61 S32-S27
      9  514  512 (S39E39 L=301)   933    125        1 ( 465/ 478- 568/ 552)  18.4 E44-E32 S40-S35
     10 1211  458 (S39W18 L=358)   311     85        0 (1173/ 428-1253/ 499)  50.2 W15-W20 S41-S35
     11 1494  450 (S37W42 L= 22)   409     55        0 (1481/ 411-1509/ 487) -59.0 W39-W45 S39-S34
  

    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