Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1436 1486 (N36W19 L=174)  3442    742       12 (1118/1336-1725/1630)  -5.6 E07-W44 N30-N39
      2  782 1529 (N24E27 L=128)   248     77        0 ( 756/1502- 814/1559)   5.5 E29-E23 N23-N25
      3 1232 1263 (N17W08 L=162)   482    171        3 (1190/1203-1308/1299)  45.0 W05-W12 N12-N19
      4 1258 1118 (N09W12 L=167)   277    108        0 (1204/1101-1315/1142)   8.1 W08-W15 N07-N09
      5  972 1040 (S02E03 L=151)   302     69        1 ( 926/1013-1006/1069)  45.0 E05-E01 S03-N00
      6  376  753 (S30E40 L=114)   323     85        0 ( 355/ 721- 404/ 779)  53.1 E42-E37 S32-S28
      7 1007  756 (S18W04 L=159)   190     36        1 ( 982/ 737-1034/ 775)   0.0 W02-W06 S18-S16
      8 1534  712 (S10W39 L=193)  1316    307        3 (1412/ 676-1634/ 765)  25.0 W30-W45 S12-S04
      9 1758  662 (S08W60 L=215)   692    108        0 (1719/ 635-1786/ 696)  51.3 W56-W63 S10-S05
     10 1416  406 (S31W43 L=198)   658    147        0 (1363/ 387-1473/ 437)   5.7 W37-W47 S32-S29
     11 1141  272 (S48W32 L=186)   328     44        1 (1120/ 253-1165/ 293)  76.0 W29-W33 S49-S45
  

    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