Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  616 1808 (N31E55 L=122)   789    118        2 ( 590/1782- 642/1837) -59.0 E58-E51 N28-N33
      2 1295 1667 (N40E02 L=176)   404     43        0 (1267/1636-1323/1693) -21.8 E04-W01 N37-N41
      3  642 1634 (N19E42 L=136)   473     92        0 ( 600/1609- 690/1652)  38.7 E44-E38 N16-N21
      4  735 1515 (N14E30 L=147)   447    130        1 ( 708/1458- 770/1574) -49.4 E34-E26 N10-N17
      5 1493  995 (N04W27 L=205)   424    125        2 (1439/ 972-1560/1012)  20.6 W23-W31 N01-N05
      6 1592  872 (N00W37 L=215)   330     86        1 (1556/ 853-1617/ 902)   0.0 W34-W39 S01-N02
      7  735  826 (S25E12 L=165)   628    190        1 ( 628/ 794- 816/ 862)  32.5 E19-E07 S28-S21
      8 1616  759 (S05W43 L=220)   770    125        3 (1589/ 694-1640/ 826) -60.3 W40-W45 S08-S01
      9 1359  682 (S17W29 L=207)  1668    398        1 (1169/ 627-1525/ 768)  36.4 W17-W36 S22-S07
     10 1563  697 (S10W42 L=219)   394     96        1 (1540/ 647-1584/ 739) -68.2 W39-W43 S12-S07
     11 1021  448 (S40W19 L=196)  1925    373        5 ( 847/ 415-1202/ 485)  15.1 W05-W31 S44-S35
     12 1544  387 (S26W57 L=234)  1814    404        2 (1391/ 357-1654/ 424)  30.1 W44-W65 S31-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