Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  626 1667 (N22E44 L=238)  1606    177        4 ( 544/1618- 721/1703)  28.6 E49-E37 N17-N25
      2 1539 1503 (N35W21 L=303)  1296    242        2 (1424/1440-1633/1557)   3.2 W11-W29 N33-N37
      3  293 1504 (N06E64 L=218)   980    117        2 ( 246/1455- 324/1546)  80.5 E68-E61 N02-N08
      4  422 1435 (N03E48 L=234)   302     47        0 ( 411/1400- 436/1462) -76.0 E49-E46 N01-N05
      5 1800 1210 (N25W47 L=329)  1176    261        3 (1769/1103-1825/1338) -81.9 W43-W50 N17-N32
      6  194 1271 (S09E65 L=217)   659     93        0 ( 176/1228- 223/1316) -32.0 E68-E59 S11-S06
      7  823 1220 (S01E16 L=266)   367     98        2 ( 770/1192- 862/1242)   9.5 E19-E13 S02-N00
      8 1725 1140 (N18W41 L=323)   278     30        1 (1703/1125-1753/1159)  11.3 W38-W43 N16-N19
      9  819  870 (S21E08 L=274)   230    107        0 ( 772/ 843- 861/ 893)  -8.1 E11-E04 S21-S19
     10  981  770 (S22W04 L=286)   294    115        1 ( 959/ 726-1011/ 827) -38.7 W01-W06 S24-S19
     11 1845  807 (N04W61 L=343)   389     73        0 (1803/ 789-1877/ 827)  26.6 W56-W64 N02-N06
     12  544  559 (S47E23 L=259)  1182    249        3 ( 442/ 510- 652/ 608)   0.0 E33-E12 S48-S43
  

    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