Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  773 1735 (N59E10 L= 50)  1213    278        4 ( 647/1697- 905/1771)  -1.8 E26-W05 N55-N60
      2 1050 1436 (N31W08 L= 68)   364    116        1 ( 993/1393-1081/1461) -56.3 W04-W10 N27-N33
      3  840 1366 (N30E07 L= 53)   355    124        0 ( 797/1337- 902/1406) -45.0 E10-E03 N26-N32
      4 1554 1238 (N11W37 L= 97)   337     60        0 (1519/1214-1585/1261) -26.6 W35-W39 N08-N12
      5  240 1198 (N25E58 L=  2)   440     91        1 ( 221/1143- 260/1234) -50.2 E60-E55 N21-N27
      6  403 1130 (N20E41 L= 19)   421    108        1 ( 389/1081- 420/1178) -80.5 E42-E39 N16-N22
      7 1913  964 (S13W69 L=129)  2452    179        4 (1883/ 891-1946/1059) -74.7 W63-W76 S18-S07
      8 1549  956 (S06W32 L= 92)   289     73        0 (1537/ 917-1562/ 998)  78.7 W31-W32 S08-S03
      9  177  951 (N10E65 L=355)   490     55        1 ( 155/ 924- 195/ 976)  31.0 E68-E62 N08-N11
     10 1210  747 (S13W08 L= 68)   332    112        0 (1170/ 724-1261/ 783)  18.4 W05-W11 S14-S11
     11 1770  632 (S32W53 L=113)  1465    203        0 (1689/ 565-1855/ 683)   7.8 W44-W66 S35-S29
     12  385  587 (S14E52 L=  8)   833    171        1 ( 341/ 507- 408/ 666)  42.5 E60-E47 S19-S08
     13  241  603 (S12E69 L=351)   536     69        1 ( 217/ 574- 264/ 619) -31.0 E72-E64 S14-S11
     14 1274  362 (S42W08 L= 68)   501     97        1 (1232/ 334-1327/ 380)   0.0 W04-W12 S43-S39
     15 1086  268 (S47E10 L= 50)   998    126        1 (1033/ 231-1155/ 292)  13.0 E15-E02 S49-S43
     16  913  203 (S50E29 L= 31)  1258    133        1 ( 876/ 158- 952/ 256)  37.6 E36-E22 S55-S44
  

    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