Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  853 1516 (N19E24 L= 26)   249     69        0 ( 821/1499- 895/1533)  11.3 E26-E21 N17-N19
      2 1413 1427 (N30W13 L= 63)   422    151        0 (1371/1371-1461/1454)  63.4 W10-W15 N25-N33
      3 1813 1024 (N18W53 L=104)  3176    409        9 (1745/ 900-1876/1194) -58.0 W46-W61 N11-N28
      4 1610 1139 (N18W33 L= 83)   245     48        0 (1584/1123-1639/1156)   0.0 W30-W35 N17-N19
      5  132 1022 (S26E72 L=338)   853     82        1 ( 112/ 974- 148/1056) -45.0 E77-E69 S28-S23
      6  246  919 (S30E54 L=357)   599    147        0 ( 199/ 880- 295/ 948)  28.6 E60-E48 S33-S27
      7 1686  607 (S08W57 L=107)   331     51        0 (1661/ 583-1708/ 630)  76.0 W55-W58 S10-S06
      8  483  553 (S49E28 L= 23)   369     41        1 ( 460/ 535- 506/ 580) -45.0 E30-E24 S50-S46
      9  820  506 (S41W03 L= 54)   255     64        0 ( 793/ 487- 853/ 533)   0.0 W00-W05 S41-S39
     10 1031  498 (S35W18 L= 68)   513     70        1 ( 992/ 479-1079/ 524)  38.7 W14-W20 S36-S32
     11 1365  343 (S33W50 L=100)   461     82        0 (1318/ 325-1394/ 362)   6.3 W44-W53 S34-S31
     12 1144  170 (S53W57 L=108)   557     57        1 (1119/ 144-1163/ 202) -22.6 W51-W64 S54-S49
  

    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