Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  905 1851 (N57E34 L= 85)   611     95        1 ( 860/1835- 956/1866)  11.3 E39-E27 N55-N58
      2  814 1691 (N41E29 L= 90)   365     66        0 ( 778/1662- 843/1710) -15.9 E32-E25 N38-N42
      3  699 1623 (N34E34 L= 84)   349     95        0 ( 671/1584- 731/1663) -32.0 E38-E30 N31-N36
      4  983 1625 (N38E13 L=106)   596    145        1 ( 931/1595-1065/1644)   0.0 E16-E06 N35-N39
      5 1472 1395 (N29W26 L=145)  2300    471        4 (1253/1306-1663/1462)  -5.2 W09-W42 N25-N31
      6  227 1112 (S04E58 L= 60)   450     92        0 ( 201/1071- 261/1163)  50.2 E60-E55 S07-N00
      7 1665 1132 (N15W42 L=160)   394    142        0 (1648/1067-1707/1163) -80.5 W40-W45 N10-N16
      8  440 1093 (S03E38 L= 80)   325     81        0 ( 422/1052- 459/1139)  90.0 E39-E37 S05-N00
      9  674 1064 (S02E22 L= 97)   321     81        0 ( 638/1039- 711/1090)  45.0 E23-E19 S03-N00
     10 1008  767 (S15W02 L=121)   230    105        0 ( 976/ 727-1024/ 805) -63.4 W00-W03 S17-S13
     11 1717  692 (S11W54 L=172)   650    142        0 (1667/ 639-1760/ 753)  63.4 W50-W56 S14-S06
     12 1679  492 (S24W61 L=179)  1478    212        2 (1564/ 468-1759/ 523)  13.4 W49-W72 S26-S21
  

    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