Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  781 1671 (N51W04 L=264)   457     69        2 ( 751/1638- 806/1709)  39.8 W01-W08 N48-N53
      2 1454 1678 (N29W52 L=313)   416     64        1 (1421/1656-1488/1694)  -9.5 W49-W55 N27-N30
      3  295 1191 (N31E47 L=213)   633     99        2 ( 250/1153- 341/1223)  11.3 E52-E42 N29-N32
      4  653 1128 (N19E19 L=241)   382    109        0 ( 608/1095- 701/1168) -50.2 E21-E16 N16-N22
      5  423 1087 (N22E36 L=224)   331     83        0 ( 386/1071- 473/1108) -23.2 E39-E32 N20-N24
      6 1372  740 (S20W13 L=273)   239    115        0 (1343/ 692-1411/ 782)  26.6 W10-W16 S22-S18
      7 1680  700 (S33W36 L=296)   285     29        0 (1658/ 684-1706/ 713)   0.0 W33-W38 S33-S31
      8 1397  643 (S27W12 L=273)   264     72        2 (1364/ 616-1421/ 671)  11.3 W09-W14 S28-S25
      9  878  383 (S29E28 L=232)   243     49        1 ( 860/ 366- 901/ 408) -76.0 E29-E26 S31-S27
     10 1445  284 (S58W09 L=270)   421     45        0 (1417/ 271-1472/ 305)   0.0 W06-W12 S60-S56
     11 1267  192 (S60E14 L=246)   634     59        0 (1243/ 171-1302/ 220) -66.8 E17-E11 S63-S55
     12 1144  157 (S58E29 L=231)  1524    166        1 (1059/ 130-1220/ 190) -17.4 E40-E20 S62-S53
  

    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