Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  660 1791 (N35E52 L= 94)   906    241        2 ( 617/1722- 700/1849) -68.2 E57-E48 N29-N40
      2 1047 1732 (N42E21 L=125)   761    126        1 (1015/1681-1076/1779) -45.0 E25-E17 N37-N45
      3  696 1606 (N23E37 L=109)   297     44        0 ( 676/1579- 716/1637) -76.0 E38-E36 N21-N25
      4  451 1547 (N14E53 L= 94)  1849    336        4 ( 364/1463- 521/1598)  55.0 E60-E45 N08-N18
      5  613 1504 (N15E38 L=108)   411     83        2 ( 590/1460- 642/1539) -68.2 E39-E36 N12-N17
      6  879 1386 (N14E17 L=129)   462    183        0 ( 834/1336- 926/1415)  18.4 E20-E14 N10-N16
      7  734 1316 (N07E24 L=122)   293     87        1 ( 713/1266- 750/1353) -80.5 E25-E22 N03-N09
      8 1117 1057 (N01W04 L=151)   288     94        2 (1091/1029-1143/1083)  56.3 W02-W05 S01-N02
      9  921  833 (S17E01 L=145)   277     56        1 ( 888/ 798- 949/ 860) -45.0 E03-W00 S18-S15
     10 1730  527 (S13W66 L=212)  1575    215        2 (1662/ 468-1802/ 596)  45.0 W61-W71 S17-S07
     11  674  397 (S53E08 L=138)   326     76        0 ( 634/ 383- 714/ 415)   7.1 E12-E04 S54-S51
     12 1125  370 (S40W28 L=174)  1024    239        1 (1021/ 342-1252/ 399)  12.5 W19-W37 S42-S37
  

    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